Decoherence of charge states in double quantum dots due to 

cotunneling 

U. HartmannQ and F.K. Wilhelm 

Sektion Physik and CeNS, Ludwig-Maximilians-Universitat, 
Theresienstr. 37, 80333 Miinchen, Germany 

Abstract 

Solid state quantum bits are a promising candidate for the realization of a scalable quantum 
computer, however, they are usually strongly limited by decoherence. We consider a double quan- 
tum dot charge qubit, whose basis states are defined by the position of an additional electron in the 
system of two laterally coupled quantum dots. The coupling of these two states can be controlled 
externally by a quantum point contact between the two dots. We discuss the decoherence through 
coupling to the electronic leads due to cotunneling processes. We focus on a simple Gedanken ex- 
periment, where the system is initially brought into a superposition and then the inter-dot coupling 
is removed nonadiabatically. We treat the system by invoking the Schrieffer- Wolff transformation 
in order to obtain a transformed Hamiltonian describing the cotunneling, and then obtain the 
dynamics of the density matrix using the Bloch-Redfield theory. As a main result, we show that 
there is energy relaxation even in the absence of inter-dot coupling. This is in contrast to what 
would be expected from the Spin-Boson model and is due to the fact that a quantum dot is coupled 
to two distinct baths. 
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FIG. 1: Sketch of the double dot system. The coupling of the double dot to the leads is assumed 
to be weak, whereas the coupling between the dots can be strong. The leads are biased such that 
sequential tunneling is supressed. 

Quantum dots ("artificial atoms") are prototype systems for studying the properties of 
discrete levels embedded in a solid-state environment [|TJ. In particular, various schemes for 
realizing quantum bits, fully controlled quantum coherent two-state systems, using quantum 
dots have been brought forward. Next to using optically excited charge states in quantum 
dots @ and electronic quantum dots used for spin manipulation H, it has been proposed |4| 
to use the charge states of a double quantum dot as a computational basis. The proposed 
setup is sketched in Figure |l|. In order to minimize the inevitable decoherence through 
coupling to the electronic leads, the system can be brought into the Coulomb blockade 
regime where sequential tunneling is supressed. We are going to discuss in this article, 
how the inevitable cotunneling still decoheres the system in this regime. The calculation 
is carried out for one specific Gedanken experiment which should capture the most generic 
features, the decay of a superposition state when the coupling between the dots is switched 
off. A more complete treatment of this setup is in preparation ||. 

We restrict our analysis on spin-polarized electrons. The relevant Hilbert space is char- 
acterized by four basis states, written as which denotes i additional electrons on the 
left dot, j additional electrons on the right dot. The two states |1,0) and |0, 1) define the 
computational basis ||. In order to describe cotunneling, we use the closest energetically 
forbidden states as virtual intermediate states. These are \vq) = |0,0) and [1)2) = 
Zero- and two electron states with internal polarization are energetically even less favorable 
due to the high charging energy of the individual dots. 
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The Hamiltonian of this system can be written as 

H = H + H 1 (1) 

H = e as (n, - h r ) - e a h vo + €ph V2 

+ 72>?<tf + «) + E + E « < 2 > 

#1 = *c + + to £ (O^frg + (3) 

Note, that the sum over dot states only runs over the restricted Hilbert space described 
above. H describes the energy spectrum of the uncoupled system, whereas the tunneling 
part Hi describes the coupling of each dot to its lead and will be treated as a perturbation. 
hi/ r are the number operators for the additional electrons on either dot. The asymmetry 
energy e as describes the difference between the energy level for the additional electron in left 
dot and the corresponding energy level in the right dot. It can be tuned through via the 
gate voltages which are applied at each dot. and e a are the energy differences towards 
the higher level j^) and the lower level \vq) respectively. 7 is the tunable inter-dot coupling. 
The a^s and b^s denote the creation/destruction operators in the dots and leads. In Hi 
the symbol t c represents the coupling constant concerning the coupling of the dots to the 
leads, which should be small compared to the asymmetry energy. Note, that we have chosen 
a slightly asymmetric notation in order to highlight the physical model: For the actual 
calculation, Hi is also expressed in the eigenstate basis of the dot. 

For our Gedanken experiment, we assume that first the inter-dot coupling 7 is large 
(7 ^ £a,s,V) such that the system relaxes into the ground state, which is a molecular 
superposition state of the form \g) = (|0, 1) — |1, 0))/y/2. Then the gate voltage that controls 
the inter-dot coupling is switched to high values, so that the coupling is practically zero. 
After this, the system dephases and relaxes into a thermal mixture of the localized eigenstates 
of the new system. 

Thus, in order to describe decoherence, we only have to consider the case 7 = K. This 
means, that H is already diagonal, i.e. the states |1, 1), |1, 0), |0, 1) and |0, 0) are eigenstates 
of our system. 

The decoherence is analyzed applying the well-established Bloch-Redfield theory, which 
is based on the Born approximation in the system-bath coupling. As we are in the Coulomb 
blockade regime, the rates evaluated from the original coupling Hamiltonian Hi vanish in 
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that order. In order to treat cotunneling with this formalism, we perform a generalized 
Schrieffer-Wolf transformation J7|, [8]]. This transformation maps our original Hamiltonian 
Hi, which is zero in the computational basis but couples the computational states to the 
lwo/2) onto a Hamiltonian which does not have this coupling to higher states but which has 
nonzero matrix elements in the computational basis. The new terms in the Hamiltonian 
describe the amplitude of transitions between the basis states via the intermediate states. 
We perform this transformation perturbatively up to second order, i.e. all processes involving 
at most one intermediate state are taken into account. 

The new Hamiltomian Hi in our special case then can be written as 



Hi,- 

Hr- 



A{R], R, ++) bX + A(L, Lb ++) 6f6f + 
A{L], L, — ) btfb? + A(R, R\, — ) bX 
A(R], L, +-) b$tf + A(L, R\, +-) bt b% 
A{L\, R, -+) ftjtftj + A{R> Lf) _ +) b R b ^\ 



(4) 
(5) 
(6) 
(7) 



The + and — signs are indizes for the states |1, 0) resp. |0, 1). We call the As Schrieffer- Wolff 
coefficients, they are calculated along the lines of using mainly second order perturbation 
theory. For example, A(R\, R, ++) is 



A{R\,R,++) = ^ 



1 



(8) 



We now use the Bloch-Redfield equations || [10 



Pnm(t*) i^nmPnm ^ ] RnmklPklit*) 

k,l 



(9) 



where R nm ki are the elements of the Redfield tensor. These equations of motion for the 
reduced density matrix are obtained within Born approximation in the effective system- 
bath coupling, so after the Schrieffer- Wolff transformation, R is of order t c . Let us remark 
that our perturbation theory naturally breaks down below the Kondo temperature, which 
can however be made arbitrarily small by lowering t c through pinching off the contacts to 
the reservoirs. 

The Bloch-Redfield equations are of Markovian form, however, by properly using the free 
time evolution of the system, they take into account all bath correlations which are relevant 
within the Born approximation [JIT] . 
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The Redfield tensor has the form 
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(10) 



The rates entering the Redfield tensor elements are given by the following Golden-Rule 
expressions 
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(11) 
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where Hi appears in the interaction representation (written as Hi). In our formalism, it is 
of crucial importance that the expectation values over Hi vanish, i.e. that the bath produces 
only noise. As a number of expectation values of Hi turns out to be finite, we tacitly replace 



Hi by Hi — {Hi) in eqs. (|TT|) and ([12]) and use the finite expectation values to renormalize 
the diagonalized, unperturbed Hamiltonian H — > H + (Hi). In our case, the effect of this 
renormalization is of the order of 0.1 % of the original matrix elements of Hq. 
After a straightforward calculation of the above Golden-Rule rates, one gets in the general 
case a (large) sum over terms with the generic form 
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where c 



O ne can express the coupling to the leads t c by t c 



^2 ■ — , where q 

S7T n ' v 



is a conductance in terms of the quantum conductance, Ep is the Fermi energy of the leads 
and n is the number of electrons in the leads. Consequently, c is then changed to c = g|^- 



The e a and e b are terms containing varying sums or differences of ep, e a and e as . Due to the 



multitude of possibilities for virtual transitions, each element of the Redfield tensor contains 
a number of terms of this generic structure. 

In the above equations, the terms containing the Fermi function /(e) only play a role 
close to resonance and can be neglected in the Coulomb blockade []E|. The n;/ r 's represent 
Bose functions for the electron-hole pairs (excitons) that are generated during the virtual 
processes. The ?/>'s denote Digamma functions and hence diverge logarithmically at low 
temperatures. 

By solving equation (^), one finds that the off-diagonal elements decay towards zero on 
a time scale (dephasing time) whereas the diagonal density matrix elements equilibrate 
on a time scale r r (relaxation time). 

Using the above expressions, we find the rates as 
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Z is a function containing several -^-functions (or logarithms). Y\, Y-\ and ^1,-1 are different 
functions of several -^'-(Trigamma-) functions (or reciprocals), however, these functions only 
have a very weak temperature dependence. The most important part of the temperature 
dependence comes in through \ and in ni/ r and is summarized in Figure |2|. We find in 
Figure ||] that the temperature dependence is similar to the Spin-Boson case |13[ . This can 
be confirmed by inspection of the formulas (|17|)-(^TD: For the relaxation rate, one has only 
Bose functions taken at the finite amount energy which is dissipated. In case of the dephasing 
rate, there are also terms that are proportional to T, which represent dephasing processes 
which do not change the energy of the qubit, i.e. cotunneling processes which originate and 
end in the same state. This explains the observed behaviour. Note, that in the Spin-Boson 
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FIG. 2: Relaxation and dephasing time (r r and t^) as a function of temperature T, with /ij = 
0.85 K, fi r = 0.91 K, e as = 0.1 K, ep = 11 K, e a = 9 K, g = 0.1, V = 10~ 12 m 2 , £ F = 5 meV and 
n/V = 1.7- 10 15 m~ 2 . 

case, where there is only one lead, the situation corresponding to our Gedanken-experiment 
(no tunneling between the classical states) would correspond to pure dephasing, whereas in 
our system relaxation is always possible by extracting an electron on one side and adding 
one on the other side from the other lead. 

The numerical values for the relaxation and dephasing times are comparedly huge, on 
the order of 100 milliseconds as compared to the experimentally measured times, which 
are in the order of nanoseconds. Other possibilities to explain the small decoherence time 
are phononic and/or photonic baths |14|, [15|, [U|, or the influence of the whole electronic 
circuitry. 

We analyzed relaxation and dephasing processes in a system of two laterally coupled 
quantum dots which is coupled to two electronic (i.e. fermionic) baths. We showed that 
even in the case of vanishing inter-dot coupling, the system's energy can relax, unlike in the 
Spin-Boson model. On top of that, the temperature dependence of the rates resembles that 
of the Spin-Boson model. We identify, that this originates in the fact that the cotunneling 
rates are mostly sensitive to the distribution function of excitons. 

As a next step, the case where the inter-dot coupling 7 has finite values will be considered 
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